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Under quite general conditions critical phenomena can be described with high 
order linked cluster expansions. The coefficients of the series admit a graphical 
expansion that is generated with the aid of computers. Our generalization of 
linked cluster expansions from an infinite to a finite volume allows to perform a 
finite size scaling analysis. We also indicate a generalization to Dynamical Linked 
Cluster Expansions with possible applications to spin glasses and neural networks 
with coupled spin and interaction dynamics. 



1 Linked Cluster Expansions in the Infinite Volume 

We will focus on multiscale phenomena that occur as critical phenomena in 
statistical systems at second order phase transitions. We consider statistical 
systems with regions of first and second order transition regions in phase space. 
Our computational tools are analytical rather than numerical calculations. We 
use convergent series expansions in a parameter called hopping parameter k. 
When applied to a calculation of the free energy and connected correlations, 
the expansions amount to Linked Cluster Expansions (LCEs). The expansion 
coefficients admit a representation as a sum of connected graphs. The price 
one has to pay for describing critical phenomena with an accurracy of the order 
of 1% is a high order in the expansion parameter. Since the number of con- 
tributing graphs rapidly increases with increasing order in k, a computer aided 
algorithmic generation of graphs becomes unavoidable. Even an optimization 
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of algorithms becomes essential to make the handling of millions of graphs 
feasible. Such a large number of graphs has to be dealt with to compute two 
and four-point functions to 20th order in k in LCEs. For example, this order 
is taken into account in scalar 0{N) models at finite temperaturejn order to 
measure the tiny finite temperature effects on the phase structure EJ. 

Originally, linked cluster expansions have been developed in the infinite 
volume. The reason why we have generalized them to a finite volume was 
twofold. One reason was to identify f st order transitions from series expansions 
in the high temperature phase. The second reason was to distinguish second 
order transitions associated with different universality classes. 



1.1 Algorithmic Generation of Graphs 

A convenient representation of any graph T is provided by an incidence matrix 
Jr with matrix elements Ir{i,j), i,j^ 1 . . .V , V being the number of vertices 
of r. Let us enumerate the vertices in some way. Ir{i,j) is then given by 

Ir{i,i) = number of external lines at (f) 
^r(i,i) = number of lines connecting vertices Vi and Vj, i,jE 1 . . .V (2) 

This representation is not unique. For any permutation of vertices, /r will 
change unless it corresponds to a symmetry of F. A possible way out of this 
ambiguity is the introduction of a complete order relation among the Is. (For 
instance, /(^^ > I^^) if {i'J') exists such that I^^\i',j') > I^'^\i',j'), and 
I^"^Hhj) — I'^^Hhj) for all i = i' , j < j' and for all i < i' and arbitrary 
J.) A unique representation of the graph is then defined as the maximum 
over all permutations of the vertices of F, i.e. as /p'*'^ = max^rePv ^r.? where 
-^f j) = -^r(7r(z), 7r(j')). For the simple example of the graph of Fig. |l| Ir and 



Figure 1: 



Example for a graph with 2 external and 6 internal lines with labels attached to 
the vertices. 
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jmax g^j.g given by 
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A CPU-time consuming part in computer aided LCEs is the complete and 
unique generation of graphs. It is performed iteratively with a certain algo- 
rithm for adding new lines and new vertices. It turns out that in order to make 
the algorithm both complete (i.e. generating all graphs that contribute) and 
efficient, multiple generations cannot be completely avoided. Let us assume 
that for a particular class of graphs we have just generated the n + 1st graph 
after the first topologically inequivalent n's have been generated. The time 
consuming procedure is then to find out whether the n + 1st graph already 
exists among the first ns or not. One has to run through all /r^ ,i = 1, . . ■ ,n 
to decide whether r„+i is topologically inequivalent to each of the first ns. A 
comparision between any 2 graphs is cheap if they are uniquely represented. If 
it is achieved via I^^^ as indicated above, V\ permutations of vertices have to 
be performed for every newly generated graph to obtain this representation. 
This is quite a large number if one goes to a high order in the expansion pa- 
rameter (for example it is 11! for the graph of order with 12 internal lines 
shown in Fig.^j.) 

To avoid the high factorials of permutations of incidence matrices, one 
needs a refined representation of graphs. One possibility is to first (partially) 
order the vertices and then to perform a maximization procedure. A vertex 
ordering that takes all local pmperties of vertices into account has been intro- 
duced by Liischer and Weisz U. With this improvement the 14th order in k 
ecame feasible. A further refinement was an extended ordering proposed in 
which in addition accounts for nonlocal topological properties to distinguish 
vertices in combination with an iterated ordering. The type of ordering of 
vertices (if there is any) is reflected in the labelling of the vertices, cf. Fig. ||. 



Figure 2: Different levels of ordering of verticesi-i a. no ordering, b. according to | 

according toEI. 

3 3 3 4 5 3 




(a) ^ 9! (b) (3!)3 (c) 

Fig. H a shows a graph F with na labelling corresponding to no ordering, 
Fig. 1^ b with^ ordering according tou and Fig. c with an extended ordering 
according to El. The advantage now is that in order to define a final unique 
representation by means of the above maximization procedure, one has to 
perform permutations only among vertices with equal labels, this is 11! for 
Fig. H a, 9! for Fig. | b, but only 3!^ for Fig. | c. At best the number of 
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remaining permutations is as small as the symmetry factor of the graph, i.e. 
the number of simultaneous permutations of rows and columns that leave Ir 
invariant. It turned out that the gain in CPU-time that is achieved because of 
the reduced number of permutations in / is considerably larger than the loose 
because of the ordering of vertices. Further detailsLon the (extended) ordering 
and the algorithmic performance can be found inEJ. 

2 The Monotony Criterion and its Application to Monte Carlo Sim- 
ulations 

o 

In LCEs the action is split into a sum of ultralocal parts S and a next neighbour 
part Snn with next neighbour couplings oc k. A Taylor expansion in n of the 
logarithm of the partition function In Z about the ultralocal contribution to 
In Z finally leads to graphical expansions of n-point susceptibilities Xn with 
coefficients a^^^ For every /i, n, aj^'^ is a sum over all connected graphs with 
/i internal and n external lines each of which adds as its weight a product of 
the inverse topological symmetry factor, an internal symmetry factor, a lattice 
embedding factor, and a product of vertex contributions depending on the 

o 

couplings involved in S- It is only the embedding factor that depends on the 
topology of the particular lattice that will change in passing from an infinite 
to a finite volume. 

For a certain interval of the scaling region response functions with a nonan- 
alytic behaviour in the infinite volume limit show different monotony behaviour 
for 1st and 2nd order transitions. Examples for such functions are the specific 
heat and order parameter susceptibilities. They are increasing in volume in a 
certain neighbourhood of Tc for 2nd order transitions, and decreasing for 1st 
order transitions for some range in the scaling region, which has to be further 
specified. For definiteness we fix the notation in terms of order parameter 
susceptibilities x- In particular, at Tc{V) , x has a "5- function" or power law 
type of singularity for a 1st or 2nd order transition in the thermodynamic limit, 
respectively. It is this difference that is responsible for the different monotony 
properties in the finite volume. Inl3 we have made these statements more pre- 
cise in order to show that the monotony behaviour of the susceptibilities is 
neither a pecularity of specific models nor an artifact of the series expansion. 
It is a generic feature of models with first and second order transitions if the 
standard assumptions on their finite size scaling behaviour apply. Here we 
further comment on applications to Monte Carlo calculations. 

Choose two volumes V^i, V2 both sufhciently large in order to satisfy the 
standard assumptions of a FSS so that the regular contribution to x and its 
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induced generic volume dependence are neglegible. Let Vi ^ V2 as argued 
below. Let A denote a generic coupling parametrizing the transition surface 
in coupling constant space. Choose S such that cia{V2)~ < \S\ < C2(j{Vi), 
where Ci , C2 and e are positive constants and a denotes the width of the critical 
region. Define 

The monotony criterion says that 

{> , 2nd order 
< , 1st order (5) 
= , tricritical point for dr/dX ^ 0. 

The difference of using the monotony criterion in Monte Carlo simulations or 
in series representations of x comes from the second volume V2 > Vi that must 
be finite in Monte Carlo simulations, but can be infinite in LCEs. For V2 < 00 
also in case of 1st order transitions a small neighbourhood around the peak of 
X exists where x increases with V because of the rounding of the (5-singularity. 
Thus V2 should be chosen sufficiently larger than Vi so that the widths of the 
critical regions behave like a{V2) ""(Vi). For finite a{V2) the points T or k 
at which the xs are evaluated should be chosen at the same distance from the 
(volume dependent) position of the peak of x- 

Deviations from these predictons in Monte Carlo calculations may be 
caused by • Vi,V2 both too small • V2/V1 too small • (k — HciV)) / i^dV) 
too large, i.e. contributions from the regular part cannot be neglected. • 
(k — Kc{V)) / Kc{V) too small, so that r > even for 1st order. 

In particular the interval of allowed ks in Eq.(^ depends on the coupling A. 
Note that the ratios of Eq. (^) do not include the generic volume dependence 
that is induced by the analytic part of x- Assuming that one is lucky in 
simultaneously matching these constraints, a result of ry^y^ > for all k 
from the specified scaling region excludes (indicates) a 1st order transition, 
respectively. 

If there is a crossover phenomenon rather than a true phase transition, r 
should vanish for a finite neighbourhood of A for sufficiently large Vi,V2. 

An application of the monotony criterion to LCEs in a scalar 0{N) theory 
with and <I>®-tcrms in 3 dimensions with iV = 1 or = 4 improvecLthe 
localizability of the tricritical line in this model by 2 orders of magnitudeH. 
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3 Dynamical Linked Cluster Expansions 

As Dynamical Linked Cluster Expansions (DLCEs) we call LCEs with dynam- 
ical next neighbour couplings. The next neighbour couplings still play the role 
of expansion parameters, but are endowed with their own dynamics. More 
precisely DLCEs amount to the following generalization of the familiar linked 
cluster expansions. 

o 1 

Let A denote the set of sites of a hypercubic lattice, A the set of unordered 
pairs of sites / = (xy), not necessarily nearest neighbours. We consider a 
statistical system characterized by the partition function 

Z{H, J-v)=N I n '^^^ • n '^^^ ■ [-^(*' ^) + H + H -^i^'] 

o 1 o 1 

xeA leA xeA leA 

(6) 

with normalization factor Af such that Z{0, 0,v) = 1 and an action S having 
the form 

U;v) = ^S .) + E ^ (^') - ^ E «x,$xC/(..)^, (7) 

o 1 o 

x&A ;eA x,y£A 

with Vxy = Vyx, Vxx = 0. <I> can be an iV-component scalar field associated with 
the lattice sites a; or a Z^r-spin, likely U can be any scalar, vector or tensor field 

o 

associated with lattice links I. The action S is split into two ultralocal parts S 
1 

and S-, depending on single sites and on single links via the fields and Ui, 

and an interaction part with coupling constants v,j.y. Note that formally the 
former coupling v^y is replaced by v^y • U^y with the dynamics of U governed 

by 5. 

A new type of graphical expansion of n-point correlation functions is then 
induced by Taylor expanding In Z(H, J; v) = W{H, J; v) about v = with 

W{H,J;v = 0) = J2w{H) + J2w{J). (8) 

o 1 

xeA leA 

To illustrate the new features, we list the terms occurring in the second deriva- 
tives in the Taylor expansion of W . Written as indices, let J(/„i) and Hi stand 
for the derivatives of W with respect to J^im) and Hi at v = 0, respectively. 
The expressions 
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suggest a graphical representation according to Fig. respectively. 



(9) 



Figure 3: Graphical representation of the terms of Eq. ( respectively. 





The dashed lines indicate the new type of local connectivity, originating 
from W 1 II lines are now connected if at least one of the following conditions 
is met. • The lines share a common vertex u„= aiJT as before. • The lines 

ic ^ 1 

are part of a multiple- line sharing a common factor 5„= 

Thus the building blocks of the graphical expansion for DLCEs are m-point 
vertices and n-multiple-lines as shown in Fig. |[ 



Figure 4: Graphical rules for DLCES. 
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Note that the graphs of the former LCEs are contained as a small subclass 

of the graphs of DLCEs. If we replace for constant U vU by v and choose 
1 ^ 1 

W (J) as Ji ■ V, we have = v (1-lines as before), but = for n > 1 

(corresponding to no multiple-lines). 

A rapid increase in the number of graphs contributing to an n-point sus- 
ceptibility sets in at low order in the expansion in k. For example the number 
of graphs contributing to X2 to order is of the order of 100, while it is by 
an order of magnitude smaller in the pure LCE of X2- 
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The generic form for the action in Eq. has interesting physical appH- 
cations to systems with coupled dynamics with (fast) spins and (slow) inter- 
actions in spin glasses and neural networks, in which both spins and couplings 
are endowed with a dynamical law and evolve in timeu. 

So far the graphs for LCEs could be generated on a SUN-workstation with 
a CPU-time of the order of days and a working^pace of the order of 100 MB. 
Work on computer aided DLCEs is in progress 
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